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Abstract. 4U 0142+61 is one of the Anomalous X-ray Pulsars exhibiting hard X-ray emission above 10 keV discovered with 
INTEGRAL. In this paper we present detailed spectral and temporal characteristics both in the hard X-ray (>10 keV) and 
soft X-ray (<10 keV) domains, obtained using data from INTEGRAL, XMM-Newton, ASCA and RXTE. Accumulating data 
collected over four years with the imager IBIS-ISGRI aboard INTEGRAL, the time-averaged total spectrum shows a power- 
£> law like shape with photon index T = 0.93 ± 0.06. 4U 0142+61 is detected up to 229 keV and the flux between 20 keV and 

229 keV is (15.01 ± 0.82) x 10 -11 ergcirr 2 s _1 , which exceeds the energy flux in the 2-10 keV band by a factor of ~2.3. 
Using simultaneously collected data with the spectrometer SPI of INTEGRAL the combined total spectrum yields the first 
evidence for a spectral break above 100 keV: Assuming for the spectral shape above 20 keV a logparabolic function the peak 
energy of 4U 0142+61 is 228*^ keV. There is no evidence for significant long-term time variability of the total emission from 
4U 0142+61. Both the total flux and the spectral index are stable within the 17% level (lcr). Pulsed emission is measured with 
ISGRI up to 160 keV. The 20-160 keV profile shows a broad double-peaked pulse with a 6.2<x detection significance. The 
total pulsed spectrum can be described with a very hard power-law shape with a photon index T = 0.40 ± 0.15 and a 20-150 
keV flux of (2.68 ± 1.34) x 10~ u ergcrrT 2 s . To perform accurate phase-resolved spectroscopy over the total X-ray window, 
we produced pulse profiles in absolute phase for INTEGRAL-ISGRI, RXTE-PCA, XMM-Newton-FN and ASCA-GIS. The 
two known pulses in all soft X-ray profiles below 10 keV are located at the same phases. Three XMM-Newton observations in 
|^>J ' 2003-2004 show statistically identical profiles. However, we find a significant profile morphology change between an ASCA- 

GIS observation in 1999 following a possible glitch of 4U 0142+61. This change can be accounted to differences in relative 
strengths and spectral shapes (0.8-10 keV) of the two pulses. The principle peak in the INTEGRAL pulse profile above 20 keV 
is located at the same phase as one of the pulses detected below 10 keV. The second pulse detected with INTEGRAL is 
slightly shifted with respect to the second peak observed in the soft X-ray band. We performed consistent phase-resolved 
spectroscopy over the total high-energy band and identify at least three genuinely different pulse components with different 
spectra. The high level of consistency between the detailed results from the four missions is indicative for a remarkable stable 
geometry underlying the emission scenario. Finally, we discuss the derived detailed characteristics of the high-energy emission 
of 4U 0142+61 in relation to three models for the non-thermal hard X-ray emission. 
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1 . Introduction down. Different origins of energy budgets proposed to explain 

the high luminosities were not consistent with all observa- 
Anomalous X-ray Pulsars (AXPs) are a special subset of young , ■ £ AV n , *■ [777 — Z — c to. ™n 

J \ . j a tional characteristics of AXPs (c.f. Woods & Thompson 200« 

isolated neutron stars with pulse periods in the range ~2- 

12 s and with high period derivatives on the order of 10~ 10 - 

10~ 12 s s . Their (anomalous) X-ray luminosities (Lx ~ 10 33 - 

10 35 ergs -1 within 2-10 keV) are orders of magnitude too 



Kaspi 2007, for a review and latest developments). However, 
for Soft Gamma-ray Repeaters (SGRs), which are another sub- 
set of young isolated X-ray pulsars with similar periods, a 

, , . , toroidal magnetic field in the neutron-star core as e nergy reser- 

high to be explained by rotational energy released due to spin ,, 17: . __ IJinmk t^, • 

b r J bJ r voir was proposed by Duncan & Thompson (1992). Their mag- 
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{B > 10 15 G). The magnetar model also explains numerous of 
other observable characteristics from AXPs like the large pe- 
riod derivatives and bursting behaviour. Most features of the 
magnetar model can be found in Duncan & Thompson ( 1992); 
Thompson & Duncan! dl993l ll995J Il996l) ; iThompson etafl 
(120021) . Observational evidence collected over the last five 
years clearly showed that AXPs and SGRs have multiply simi- 
lar characteris tics and therefore are now b elieved to be subsets 
of magnetars dWoods & Thom pson 2006). 

In this paper we will present new results on 
AXP 4U 0142+61. This particular AXP has a spin 
period of P = 8.69 s and a period derivative of 
P - 0.20 x 10~ n ss , The inferred surface magnetic 
field strength is B = 3.2 x 10 19 VPP G = 1.3 x 10 14 G. In 
the 2-10 keV band this is the brightest persistent AXP with a 
lumi nosity Lx ~ 1 X 10 35 ergs" 1 a ssuming a distance of 3.6 



kpc dDurant & van Kerkwiikll2006al) . 

A remarkable discovery with the hard X-ray/soft gamma- 
ray telescope INTEGRAL dWinkler et al.l2003l) is the existence 
of a hard spectral comp onent in the sp e ctra o f AXPs at energies 
above 20 keV. In 2004 lMolkovet alj d2004l) reported hard X- 
ray emission (>20 keV) from the position of AXP IE 1841-045 
in the supernova remnant G27.4+0.0 (Kes 73). Kui per et al 



(2004) showed unambiguously that the hard X-rays indeed 
originate from the AXP and not from the supernova remnant 
by detecting pulsations in the hard X-ray emission from this 
source using Rossi X-ray Timing Explorer (RXTE) PCA and 
HEXTE data. 

For two other AXPs I NTEG RAL detections have been re- 
ported. iRevnivtsev et all d2004 and Iden Hartog et alJ d2004l) 
detected 1RXS J 170849-4009 10 (hereafter 1RXS J1708-40) 
and 4U 0142+61, resp ectively. Als o for th ese and for a fourth 



AXP, IE 2259+586, Kuiper et all d2006l) discovered pulsed 



hard X-ray emission using RXTE PCA and HEXTE data. 

The INTEGRAL total-emission spectra available for 
IE 1841-045, 1RXS J1708-40 and 4U 0142+61 are all power- 
law like and hard. The spectra of these three AXPs have photon 
indiceflof 132 + 11, 1.44 + 0.45 and 1.05 ±0.11, respect ively 
dKuiper et al.l2 006). Furthermore Jden Hartog et alj d2007l) per- 
formed a multi-frequency campaign in which 4U 0142+61 was 
observed from radio to hard X-rays almost simultaneously. A 
dedicated INTEGRAL observation of 1 Ms was part of this 
campaign. This observation yielded a power-law like spectrum 
up to 229 keV with a photon index of 0.79 + 0.10 with a lu- 
minosity of 2.6 x 10 35 ergs , which is roughly 2000 times the 
maximum luminosity possible by energy release due to spin 
down. 

Interestingly, 4U 0142+61 exhibited some bursting activi- 
ties which were discovered during the RXTE monit oring of this 



source. A first burst was found on April 6, 2006 (Kaspi et al. 



20061) . Soon ther eafter, on June 25, four more bursts were de 

tected dDib et alJ 12006). During these last event s, the source 

seemed to be in an active state. Furthermore, iGavriil et al 



d2007l) discovered a very large X-ray burst on February 7, 
2007. This burst was extraordinary as it showed characteris- 



tics similar to a SGR giant flare, only less energetic. It started 
with a bright and short spike with a peak flux of ~2.5 x 10 -8 
erg cm^s -1 (2-60 keV) and in the first part of the decay- 
ing tail single pulsations of the AXP were visible. The ex- 
tended tail lasted up to half an hour. The fluence (Tgo) of this 
burst was -14.5 x 10~ 8 erg cm^s -1 (2-60 keV), which is 6-88 
times higher than the previous five bursts from 4U 0142+61. 
Moreover, the spe ctrum shows significant spectral features 
dGavriil et al.ll2008l) . 

Theoretical attempts to find physical explanations for the 
hard X-ray emission have started, but so far none of the models 
can explain all of the published properties and for some as- 
pects the data cannot discriminate between the proposed mod- 
els. More detailed observational results are required to further 
constrain the geometries of the production sites as well as al- 
lowed production processes. In the discussion (see, Sect.|5]l we 
summarise the different theoretical approaches, the promising 
progress, but also their shortcomings. 

In this work detailed hard X-ray characteristics are pre- 
sented using all available INTEGRAL data for 4U 0142+61 
to date. These include total spectra (pulsed + DC emission) 
both time-averaged (whole mission) as well as in 5 shorter 
time stretches to look for possible long-term time variability. 
We also analysed two Target of Opportunity observations per- 
formed shortly after the above mentioned burst activities to 
look for variability of any kind. Furthermore, the time-averaged 
total-pulsed spectrum as well as time-averaged phase-resolved 
pulsed spectra are presented. 

In order to study possible relations between the hard 
(>10 keV) and soft (<10 keV) X-ray bands we have 
(re)analysed archival XMM-Newton, RXTE-PCA and ASCA- 
GIS data applying parameters and selections consistent with 
the INTEGRAL analyses. 

2. Observations and analysis 

2.1. INTEGRAL 

The hard X-ray/soft gamma-ray (3 keV - 8 MeV) mis- 
sion INTErnational Gamma-Ray Astrophysics Laboratory 
INTEGRAL dWinkler et al.ll2003l) . has been operational since 
October 2002. For this work both main instruments IBIS 
and SPI are used. IBIS, Ima ger on Board the INTEGRAL 



Satellite (Ub ertini et al.l I2003I) . is a coded-mask instrument 



1 Defining the power law as F = FqXE r with Fq the normalisation 
in units of ph cirr 2 s _1 keV -1 and T the photon index 



with a low-energy detector called IN TEGRAL Soft Gamma- 
Ray Imager (ISGRI; iLebrun et all2003l) which is sensitive be- 
tween ~20 keV and ~300 keV and has a wide field of view 
(FOV) of 29° x 29° (full-width zero response). The angular 
resolution is about 12 '. The Spectrometer for INTEGRAL (SPI; 
Vedrenne et al.ll2003l) is also a coded-mask instrument which is 
sensitive between -20 keV and ~8 MeV. The FOV is 35° x 35° 
(full-width zero response). This instrument is optimised for 
spectroscopy and has a modest angular resolution of ~2°5. SPI 
is of particular use in this work for its better sensitivity to IBIS 
above -200 keV 

In its default opera tion modes, INTEGRAL observes the 
sky in dither patterns (iJensen et al.1 120031) . Typical pointings 
(Science Windows, ScWs) have integration times of ~1800- 
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Table 1. Summary of the INTEGRAL observations of 4U 0142+61. The following information is given: the revolution intervals; 
the observations time spans both in MJD as in calendar dates; the number of Science Windows; exposure times (f e xp); and the 
effective on-source exposure (£ e ff). 



Rev. 






Time span 






ScWs 


4xp (A s ) 


t eS (ks) 


142-468 


52985 - 


- 53960 


Dec. 12, 2003 - 


-Aug. 13, 


2006 


1617 


4494.4 


2365.3 



142-189 


52985 - 


-53125 


Dec. 


12, 


2003- 


- Apr. 


30, 


2004 


404 


963.2 


342.0 


202-269 


53165- 


- 53367 


Jun. 


9, 


2004- 


-Dec. 


28. 


2004 


410 


976.5 


501.6 


331-336 


53550 - 


-53568 


Jun. 


29. 


2005- 


-Jul. 


17, 


2005 


265 


867.6 


827.9 


384-396 


53710- 


- 53745 


Dec. 


6, 


2005- 


-Jan. 


10. 


2006 


435 


1334.3 


420.6 


452-468 


53912- 


- 53960 


Jun. 


26, 


2006- 


-Aug 


13, 


2006 


103 


353.8 


273.2 



454 ToO 53918-53920 Jul. 2, 2006 -Jul. 4, 2006 59 213.7 209.7 

528 ToO 54139-54142 Feb. 8, 2007 -Feb. 11, 2007 54 198.6 195.5 



3600 s. Table [TJ lists the total exposure time for the selected 
ScWs for which 4U 0142+61 was within an angle of 14?5 
from the pointing direction of IBIS, as well as the effective 
on-source exposure, which is reduced due to off-axis viewing 
angles. Due to the visibility windows of this part of the Galaxy 
for INTEGRAL operations, the observations of the prime tar- 
gets form groups of consequentive orbital revolutions (Revs). 

There are three dedicated observations on 4U 0142+61. 
The first one was performed mid 2005 and lasted for 1 Ms 
and two were shorter Target of Opportunity (ToO) observa- 
tions which were triggered by the bursting activities detected 
with RXTE on June 25, 2006 and February 7, 2007. 

The data were screened for Solar flares and erratic count- 
rates due to the passages through the Earth's radiation belts. 
After screening 1671 ScWs, resulting in 4.7 Ms total exposure 
time, were available for the analyses. 

2.1 .1 . INTEGRAL spectra from spatial analysis 

For spectral analysis the IBIS-ISGRI data of all 1671 ScWs 
are reduced using the Off-lin e Scientific Analysis ( OSA) soft- 
ware package version 5.1 (see lGoldwurm et al.l2003l) . For each 
ScW deconvolved sky images are created in 20 energy intervals 
using exponential binning between 20 keV and 300 keV. For 
every energy band time averaged count rates are determined 
by averaging the count rates, weighted by the variances, from 
all deconvolved maps at the position of 4U 0142+61, for each 
set of consecutive observations and the total data set. To con- 
vert the count rates to flux values we have used the procedure 
described by iKuiper et al.l (12006;). The measured AXP count 
rates are normalized to the Crab count rates, determined from 
Crab observations performed during INTEGRAL Rev. 102. 
The AXP fluxes in Crab units are converted into photon fluxes 
using a curved power-law shape; 



F y = 1.5703(14)x(£ y /0.06335) 



-2.097(2)-0.0082(16)xln(£ r /0.06335) 



(1) 



2 in 

,2 



Kuiper et alJ 120061) . where F 7 is expressed in 



(Eq. 

ph/(cm^ s Me V) and E y in Me V. Finally the spectra are im- 
ported in XSPEC version 12.3 (lArnaudlll996l) . which is used 
for the spectral fitting. All errors quoted in this paper are lcr 
errors, unless stated otherwise. 



The spectral analysis of SPI data is based on all (screened) 
data up to Rev. 336. The data were binned in 9 energy bins cov- 
ering the 20 ke V - 1 MeV energy range. For each of the energy 
bins, source fluxes have then been extracted using a maximum 
likelihood fitting proced ure that considers the instrumental re- 
sponse of the telescope (Knodlseder 2004). The measured flux 
is then normalized to the Crab spectrum that has been obtained 
by fitting the SPI data of revolutions 43-45. Multiplication 
with the above-mentioned curved power-law Crab model re- 
sults then into photon fluxes. 



2.1.2. INTEGRAL timing analysis 

For the INTEGRAL ti ming analys i s, we h ave followed the pro- 
cedure described by IKuiper et al.l (I2006I) . We only considered 
ISGRI data for this purpose. Data up to Rev. 396 are used 
(see TableHJ. Events are selected from non-noisy detector pix- 
els that were/could be illuminated by the source through an 
open mask element with a greater illumination factor than 25%. 
Instrumental, on-board processing and grou nd-station time de- 
lays are corrected for (Walter et al.l I2003). The event arrival 
times at the space craft are converted to arrival times at the 
Solar-system barycenter using the JPL DE200 Solar-system 
ephemeris and th e source position measured with Chandra 
(ljuett et al.ll2002l) . The barycentered arrival times are finally 
folded on an appropriate phase-connected ephemeris (see Table 
|2} created using contemporaneous RXTE data (see Sect. I2.21 i. 
The (TDB) time to pulse phase conversion taking into account 
consistent phase alignment for each ephemeris is provided by 
the following formula: 



o(o = v-(f-f Epoch )+-v-(f-f Epoch r +-v-(f-f Epoch r-o .(2) 

1 6 

The detection significances of the derived pulse profiles are es- 
timated applying the Zj; test ( Buccheri et al.lll983l) . 

In order to derive the spectrum of the pulsed emission, 
first the number of excess counts in the obtained pulse profiles 
above flat background/DC levels have to be determined. For 
this purpose, the pulse profiles for the selected energy bands 
are fitted with truncated Fourier series using two harmonics. 
The minimum of this fit is defined as the DC level. The excess 
counts are counted above this level and converted into flux units 
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Table 2. Phase-coherent ephemerides derived from RXTE-PCA monitoring data, valid for the analysed INTEGRAL observa- 
tions. 



AXP 


Start 


End 


INTEGRAL 


'Epoch 


V 


V' 


V 


O 




[MJD] 


[MJD] 


range (Revs) 


[MJD, TDB] 


[Hz] 


xl0- 13 [Hzs-'] 


xlO- 22 [Hzs~ 2 ] 




4U 0142+61 


52549 


53169 


001-203 


52726 


0.1150945693297 


-0.267038 


0.220 


0.9442 


4U 0142+61 


53251 


53619 


230-353 


53420 


0.1150929854501 


-0.263899 


0.307 


0.8906 


4U 0142+61 


53562 


53745 


334-396 


53650 


0.1150924611500 


-0.267430 


-0.319 


0.2817 



applying the same procedure for the Crab pulsar a nd normalis- 2.3. XMM-Newton 



ingto the known Crab-pulsar spectrum (Eq. 3 in iKuiper et al 

2006); 



F y = 0.4693(21)x(£ r /0.04844) 



-1.955(7)-0.0710(78)xln(£ r /0.04844) 



(3) 



Phase-resolved spectra are similarly created using only excess 
counts within selected phase intervals (see Sect. 13.3b . 



Table 3. XMM-Newton observations of 4U 0142+61. The ob- 
servations are designated A, B and C and referenced accord- 
ingly in the text. Further, the observation dates, exposures and 
Good Time Interval (GTI) exposures for MOS and PN are 
given. 



2.2. RXTE 

The Rossi X-ray Timing Explorer (RXTE) has been operational 
since 1996. For this work we use data from the Proportional 
Counter Arrays (PCA) onboard RXTE djahoda et all Il996h . 
which is a non-imaging instrument sensitive in the 2-60 keV 
energy range. The FOV is approximately 1° (FWHM). 

Regular monitoring observation s of 4U 0142+61 with the 
RXTE-PCA (iGavriil & Kaspill2002h are the basis of the timing 
results in this paper. These observations allow us to create accu- 
rate phase-connected timing solutions. We have generated three 
new phase-connected ephemerides with consistent pulse align- 
ment to cover the whole INTEGRAL-time span. The details of 
the ephemerides can be found in Table T hey are in agree- 
ment with the solution presented bv lDib et al] d2007l) . who also 
present a phase-coherent timing solution over the period 1996- 
2006. 

We hav e used the s a me R XTE-PCA data as listed in 
Table 1 of IKuiper et al.l (|2006) which were taken between 
March 1996 and September 2003 before the INTEGRAL ob- 
servations, to perform spectral timing ana lyses analogous to the 
procedures outlined by Kuip er et al.l(|20 06). We note that in the 
folding process of PCA data the application of ephemeris entry 
one of Table |2]results in pulse profiles which are shifted 0.543 



# ObsID Date Exp. GTI-MOS GTI-PN 

(ks) (ks) (ks) 



in pha se with respect to the profiles shown in IKuiper et al. 

(2006) for 4U 0142+61. 



Irrespective of the energy band, the pulse-profiles can be 
described sufficiently accurate by truncated Fourier series with 
3 harmonics above constant background levels. The excess 
(pulsed) counts above the background levels are converted into 
flux units using properly PCU expo sure weighted ener gy re- 
sponse matrices (see e.g. Sect. 3.2 of IKuiper et al]|2006l) . 

We have extracted a time-averaged total-pulsed spectrum 
and phase-resolved spectra adopting the same phase intervals 
as selected for INTEGRAL (see Sect l3.31 >. To derive unab- 
solved (soft) X-ray spectra a Galactic Hydrogen absorption 
column of Nn = 0.57 x 10 22 ctrT 2 has been used (as fixed 
value) as determined from spectral analyses of XMM-Newton 
observations (see Sect. 13. 1.4t . 



A 0112781101 
B 0206670101 
C 0206670201 



2003-01-24 
2004-03-01 
2004-07-25 



6.4 

46.5 

23.9 



5.4 

36.4 

23.0 



3.7 

27.3 

20.6 



XMM-Newton (IJansen et al.l 1200 lb has been operational 
since 2000. Onboard are three CCD cameras for X-ray imag- 
ing, namely two EPIC (European Photon Imaging Camera) 
MOS cameras ( Turner et al.l 120011) and one EPIC-PN camera 
JStriider et alJl200ll) . All cameras have a FOV of -30' and are 
sensitive in the energy range ~0.3-12 keV. 

We have (re)analysed the four publicly available data sets 
acquired by XMM-Newton between 2002 and 2004. The data 
set taken on February 13, 2002 (Obs Id 0112780301) did not 
contain usable data for this work. The summary of the re- 
maining 3 data sets can be found in Table [3] We designate 
these observations A, B and C, respectively. In this work we 
present only the results using the PN. We have also analysed the 
data from both MOS instruments and used this for consistency 
checks. The PN was in small-window mode during observation 
A and in timing mode du ring obs ervations B and C. In parallel 
studies [Rea et al] ( I2007ab and lGonzalez et al. (2008) als o anal- 
ysed these XMM-Newton observations. lGonzalez et al] (12008) 
also included private XMM-Newton data taken in 2006 and 
2007 together with Chandra and Swift data to study long-term 
variations below 10 keV. 

The three data sets in Table[3]have been analysed using S AS 
v. 7.0 and the latest calibration files that were available (Jan. 
2007). The data have been checked for solar (soft proton) flares 
by creating a light curve for each observation for events with 
energies larger than 10 keV. From these light curves the count- 
rate distribution has been determined and subsequently fitted 
with a Gaussian. Good Time Intervals (GTIs) have been created 
by allowing only time stamps from time intervals for which 
the count rate was lower than the fitted mean count rate plus 
three times the width of the distribution (3<x). For observations 
A, B and C these cutoff values were 0.082, 0.243 and 0.358 



P.R. den Hartog et al.: Detailed high-energy characteristics of AXP 4U 0142+61 



CN 

d 



cu d 



m 
o 



E 
o 

15 O 

ci d 



q 
d 




J 





J- 



T 



20 50 100 200 20 50 100 200 20 50 100 200 20 50 100 200 20 50 100 200 

Energy [keV] 

Fig. 1. Total-flux spectra of 4U0142+61 taken during INTEGRAL revolutions (from left to right); 142-189, 202-269, 331-336, 
384-396 and 452-468. Data points are binned to 4<x if the single bins are of less significance. The upper limits represent 2<x 
confidence levels. 



counts per second. PN single and double events, i.e. patterns 
less than and equal four, and an energy range of 0.3 to 12.0 keV 
have been selected. The arrival times of the selected events are 
barycentered adop ting the m ost accurate celestial position of 
4U 0142+61 dJuettetal.1120021) . 

There was no need to correct for pile up in observation A 
as the count rate from 4U 0142+61 was sufficiently low for 
the small-window mode. A circular extraction region centered 
at 4U 0142+61 with a radius of 40" was used. For the back- 
ground extraction two circular source-free regions centered at 
different locations in the small-window were chosen with radii 
of 35" and 47'.'5, respectively. 

Observations B and C were taken in timing mode. For these 
observations we chose extraction regions with a total width of 
27 pixels centered on the source, which corresponds to 11077 
close to the 95% encircled energy radius of the source. While 
the sensitive area (ARF) generator corrects for the remaining 
~5% of missed energy outside the source extraction region, it 
does not compensate for the subtraction of any source contri- 
bution in the chosen background extraction strip. We took a 
strip of 11 pixels wide starting 102'.'5 from the point-spread 
function (PSF) core. The source contribution from this part of 
the PSF can be neglected. We fitted a Lorentzian plus constant 
background to the observed one-dimensional source profile. 
The source contribution in the 1 1 pixel wide background strip 
turned out to be a mere 0.48% of the normalized PSF, which 
means taking into account the source-to-background scaling of 
27/1 1 an ignored correction of about 1.2%. 

The extracted spectra have been binned, oversampling the 
energy resolution by a factor of three and then occasionally 
rebinned to ensure a minimum of 25 counts per bin. 

For the XMM-Newton timing analysis, the events selected 
as described above are used. These events are folded using the 
appropriate ephemeris in Table [2] For the pulsed spectra the 
pulsed excess counts are extracted using the same procedure as 
described in Sect. 12. 1.21 for the INTEGRAL pulsed spectrum, 



fitting the pulse profiles with truncated Fourier series with the 
first three harmonics above a constant background. 

2.4. ASCA 



ASCA (ITanaka et alJ 119941) was operational for eight years 
from 1993. Onboard ASCA was the Gas Imaging Spectrometer 
fGIS. IOhashi et alll996l:lMakishima et alJl996l) . The GIS was 
sensitive between ~0.7-12 keV and had a FOV of ~50'. 

We hav e revis ited the spectral timing analysis by 
Kuiper et all J2006I) of the ASCA-GIS observation of 
4U 0142+61 taken in July/August 1999 to make it consistent 
with this work, i.e., in the current spectral analysis we have 
changed and fixed the Galactic absorption column to the value 
derived from the XMM-Newton data (see Sect. 13.1.4b . In the 
timing analysis the pulse phase has been redefined to be con- 
sistent with that of INTEGRAL, RXTE and XMM-Newton, al- 
lowing the same phase selections for the phase-resolved spec- 
troscopy. 

3. Results 

In this section, we will first present the INTEGRAL and XMM- 
Newton total spectra and the search for long-term variability. 
Next, we will present pulse profiles from our timing analyses 
of INTEGRAL, XMM-Newton, RXTE and ASCA. Finally, for 
the same instruments the total pulsed and phase-resolved spec- 
tra are derived consistently, revealing that genuinely different 
components contribute to the total pulsed emission. 

3.1. Total spectrum and long-term variability 

Previ ous I NTEGRAL s tudies by Iden Hartog etail d2006 , 
2007b and iKuiper et alJ (1200 6) show that hard X-rays from 
4U 0142+61 have been detected at different epochs at simi- 
lar flux levels, suggesting a stable and persistent hard X-ray 
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Fig. 2. Error contours (lcr, 2<x and 3<x) of observations 142- 
189 (in grey) and 202-269 (in black) with respect to the time- 
averaged measurement (star). 



Table 4. Summary of the single power-law fits for the 
INTEGRAL-ISGRI total-spectra of 4U 0142+61 at different 
epochs for the observations summarized in Table [TJ For each 
epoch the photon index, 20-150 keV flux and;f 2 are given. For 
the the spectrum taken in Rev. 142^-68 also the 20-229 keV 
flux is given. 



Rev. 


r 


F 20 _i5o x 10- 11 
[ergcirr 2 s~'] 


X'r (dof) 


142-189 


0.85 + 0.14 


11.34 + 0.78 


0.69(15) 


202-269 


1.21 + 0.16 


7.80+0.83 


0.66 (13) 


331-336 


0.79 + 0.10 


9.80+0.57 


1.12(16) 


384-396 


0.89 + 0.17 


9.46 ± 0.91 


0.91 (15) 


452-468 


0.96 ± 0.25 


7.8+ 1.2 


0.82 (12) 


454 ToO 


1.06 ± 0.30 


7.8+ 1.4 


0.98 (12) 


528 ToO 


0.88 ± 0.34 


9.3+ 1.8 


1.63(11) 


142-468 


0.93 ± 0.06 


9.09 ± 0.35+ 


0.94 (16) 


F 20-229 = 


(15.01+0.82) 


x 10~ u ergcirf 


V ] 



spectrum. Firstly, we will verify the stable nature of the high- 
energy emission. Then, we will present the time-averaged to- 
tal spectrum using contemporaneous IBIS-ISGRI and SPI data. 
Finally, the total spectra of XMM-Newton are derived and 
an example is shown together with the INTEGRAL time- 
averaged total spectrum, rendering a broad-band X-ray view 
of the total spectrum. 

3.1 .1 . INTEGRAL ISGRI persistent spectra and 
long-term variability 

The ISGRI data were grouped in five sets as indicated in 
Table [TJ and analysed as described in Sect. I2.1.1I The statis- 
tical quality is not equal for every spectrum, but 4U 0142+61 
is detected significantly in each group and spectra can be ex- 
tracted from each data set. The resulting five total spectra are 
displayed in Fig. [TJ 

To fit the spectra we use single power-law models. As in 
previous work they describe the spectral shapes at a satisfactory 
level. The fit results for these spectra are listed in the upper 
part of Table H] The photon indices F range from 0.79 ± 0.10 
to 1.21 + 0.16 and the fluxes in the 20-150 keV energy band 
(^20- i so ) range from 7.8 x 10~ n to 11.3 x 10~ n ergcm" 2 s _1 . 

The errors in the indices F and fluxes F20-150 are not inde- 
pendent. In order to check for possible long-term time variabil- 
ity we therefore compared the error contours of the different fits 
with respect to the error contours of the fit to the time-averaged 
total spectrum (see below in Sect. 13.1.21 . The spectra extracted 
for Revs. 142-189 and 202-269 are the most deviating with 
respect to the time-averaged spectrum. In Fig.|2]the error con- 
tours for these spectra are shown together with the best-fit val- 
ues for the time-averaged total spectrum (Revs. 142^468). The 
two spectra deviate from the average around the ~2cr level. 

The standard deviations (s = ~\Jl/(n - 1) Yii( x i ~ *) 2 ) f° r 
the power-law indices and the 20-150 keV fluxes relative to 
the weighted means are 0.16 and 1.5 x 10~ n ergcrrT 2 s , re- 



spectively. Therefore, both the power-law shape and the 20- 
150 keV flux of 4U 0142+61 are stable within 17% (lcr). 



3.1 .2. INTEGRAL ISGRI & SPI total spectra 

Using all available (Revs. 142-468, see Table □} 2.37 Ms net 
on-source exposure the time-averaged total spectrum (Fig. O 
derived from ISGRI data can be fitted satisfactorily with a sin- 
gle power-law model. The best fit parameters give a photon 
index F = 0.93 + 0.06 and a 20-150 keV flux of (9.09 ± 0.35) x 
10~ n ergcm^s -1 . 4U 0142+61 has been detected up to 229 
keV with a 4.3<x significance in the 152-229 keV band. 

Even though the upper limits at MeV energies obtained 
with the Compton te lescope COMPTEL onboard CGRO 
dden Hartog et al.ll2.006l also shown in Fig. [3]) strongly suggest 
that the spectrum has to break at a few hundred keV, the ISGRI 
spectrum shows no indication for a spectral break. The quality 
of the fit ix 2 r = 0.94; dof =16) assuming a single power-law 
model over the 20-300 keV range is good and the use of a more 
complex model with additional parameters is statistically not 
justified. 

SPI is more sensitive than ISGRI at energies higher than 
~200 keV. Therefore, we have also analysed all available SPI 
data with 4U 0142+61 in the FOV. We have extracted the SPI 
spectrum and detected 4U 0142+61 significantly above 40 keV 
up to 140 keV. The maximum detection significance is 5<x be- 
tween 80 keV and 140 keV. In the next energy bin (140-220 
keV) a marginal detection of 2.6<x is achieved. The six flux 
values between 20 and 220 keV can be fitted (y 2 - 1 .67; dof = 
4) with a power-law model with photon index F = 0.86 ± 0.20 
and a 20-150 keV flux of (8.8±1.2)xl0" n erg cirT 2 s"\ which 



is consistent with the ISGRI total spectrum (see Fig.O. 

Fitting the SPI spectrum with a single power law using all 9 
spectral bins between 20 and 1000 keV we find a poor and un- 
satisfactory fit with a^f 2 . = 2.57 for 7 dof (x 1 - 18). Therefore, 
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1 10 100 1000 1CT 

Energy [keV] 

Fig. 3. Unabsorbed total spectra of 4U 0142+61 as measured 
with different instruments; XMM-Newton (Obs B; 0.55-11.5 
keV) in black; INTEGRAL-ISGRI (Revs. 142-468; 20-300 
keV) in black with open square symbols; INTEGRAL SPI (20- 
1000 keV) in red; and CGRO COMPTEL (0.75-30 MeV) lim- 
its in black. Shown are also the best single power-law model fit 
for the ISGRI data-points (in grey) and the 'three logparabola' 
fit for the whole band (in blue). See Sect. l3.1l for details. 



we followed lKuiper et al] ( 120061) . by using a logparabolic func- 
tion; 



F = F n X 



Eq) 



(4) 



with Eq (in units keV) the pivot energy introduced to minimize 
correlations between the parameters, and Fq the flux (in units 
phcm _2 s _1 keV -1 ) at Eq. Note that this function is a simple 
power law if the curvature parameter ft is equal to zero. 

The fit of the SPI spectrum improves with A%1 - 15.77, 
which corresponds to a ~ Act improvement for one additional 
fit parameter. 

To exploit all available spectral information we fitted the 
ISGRI and SPI spectra simultaneously. Starting from a single 
power-law model we obtained a fit with a;f 2 = 1.58 for 27 dof 
(X 1 - 42.8) which can be improved significantly. The model 
described by this power-law is too high in the energy range 
with the SPI upper limits i.e. above 220 keV. Using again a log- 
parabolic function an optimum fit with a% 2 of 20. 1 is achieved. 
A A^j of 22.7 for one additional parameter translates in a 4.8<x 
fit improvement. This is the first clear detection of the spectral 
break in the total spectrum of 4U 0142+61 above 20 keV using 
contemporaneous high-energy data. The best fit parameters are 
a = 1.26 ± 0.09,/? = 0.41 ± 0.09 and F = (5.1 + 0.3) x 10~ 6 
phcm _2 s _1 keV -1 at Eq = 91.527 keV. The peak energy (for 
a vF v representation, eq. to E 2 F) corresponding to these pa- 

e 
1 (fitting three free 



rameters is E peak = 228^ ( 3 )6 keV, while the 20-150 keV flux 



Including the non-contemporaneous COMPTEL (0.75-30 
MeV) flux measurements (upper limits in Fig. [3]) in a similar fit 
yielded best fit parameters with somewhat smaller uncertain- 
ties. We found as best fit parameters: a - 1.484 ± 0.057, ft - 
0.351 +0.044 and F = (3.01 ±0.16) x lO^phcnr 2 s -i keV" 1 
at Eq = 133.718 keV. The peak energy (E 2 F) corresponding 
to these parameters is £peak = 279*$j keV. We note that this 
result is based on an assumed "logparabolic" spectral shape. 
We adopted this spectral form in order to comply with the SPI 
and COMPTEL limits. However, we do not claim that we have 
measured the spectral shape above ~200 keV in detail. In par- 
ticular, the extension of the spectrum towards MeV energies 
is uncertain. Significant measurements in the energy range be- 
tween 200 keV and 1 MeV are required to unravel the true spec- 
tral shape in the vicinity of the break energy. 

3.1 .3. INTEGRAL post-burst spectra 

Of the five AXPs for which bursts have been detected, 
4U 0142+61 is so far the only one for which a hard X-ray 
spectrum above ~ 10 keV has been observed. After two remark- 
ably different bursting events detected with RXTE-PCA we re- 
quested two INTEGRAL ToO observations in order to study 
possible correlations between the soft- and hard X-ray emis- 
sion. For both occasions ~200 ks ToO time was granted. The 
source was observed in Revs. 454 and 528 (see Table[T|i. 

We have detected 4U 0142+61 with detection significances 
in the 20-150 keV energy band of 5.4cr and 5.9cr for Rev-454 
and Rev-528, respectively. Both spectra have been fitted with 
a power-law model which described the spectra well. The fit 
parameters for the first ToO are F = 1.06 ± 0.30 and F20-150 = 
(7.8 + 1.4) x 10~ n erg cirr 2 s _1 and for the 2nd ToO the fit pa- 
rameters are T = 0.88 + 0.34 and F20-150 = (9.3 ± 1.8) x 10~ n 
erg crrT 2 s . The values are fully consistent with those for the 
time -averaged spectrurq 2 ] (see also Table |4j. 

Gonzalez et al.l (12008b reported a flux increase of (15±3)% 



in the 2-10 keV band observed with XMM-Newton coinciding 
with the burst activity in 2006-2007. It should be noted, how- 
ever, that in order to measure with INTEGRAL a significant 
(>3cr) change in the hard X-ray spectrum within a 200 ks ob- 
servation, the flux level should have changed by at least ~60%. 
Correlated changes in the hard and soft X-ray fluxes of the scale 
reported bv lGonzalez et al.l (2008) could therefore not be mea- 
sured. 

3.1.4. XMM-Newton EPIC-PN total emission spectra 

We derived the XMM-Newton EPIC-PN total (= pulsed plus 
DC) spectra as described in Sect. l2.3l (see e.g. Fig[3]for an ex- 
ample). We first investigated the absorption column (Nh) using 
the XMM observations with the best statistics (B and C) and 
concentrated ourselves on energies below ~5 keV i.e. the soft 
part of the X-ray spectrum. 



amounts: (8.97 : 
parameters). 



0.86) x 10" 11 ergcirr 2 s 



2 Note that the 1st ToO is included in the time-averaged total spec- 
trum, but we have tested the ToO with the time-averaged total spec- 
trum of the first four data sets from which it was not statistically dif- 
ferent. 
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Table 5. XMM-Newton EPIC-PN spectral fits of the total-emission spectrum of 4U 0142+61. Three logparabolae (Eq.[4ji were 
used for an empirical fit including three fixed parameters from the ISGRI/SPI 20-1000 keV spectral fit. For each observation (A, 
B and C, see Table [3j six optimum parameter values are shown, three for each logparabolic model (indicated with a subscript 
number), for the best fit (0.55-10 keV). The last row lists the unabsorbed 2-10 keV energy flux. 





A 


B 


C 


^0,1 


(6.85 ± 0.05) x 10" 2 


(5.50 + 0.02) x 10" 2 


(4.57 ± 0.02) x 10" 2 


u-i 


0.933 ± 0.012 


0.673 ± 0.006 


0.528 ± 0.008 


A 


1.534 + 0.013 


1.899 + 0.007 


2.089 + 0.010 


■F(),2 


(1.04 + 0.04) x 10" 2 


(2.33 ± 0.02) x 10" 2 


(3.20 + 0.02) x 10" 2 


a% 


2.095 ± 0.041 


1.671 + 0.009 


1.440 + 0.008 


h 


0.212 + 0.025 


0.622 ± 0.007 


0.808 ± 0.006 


^2-iokeV (ergcnrV) 


(6.68 + 0.40) x 10-" 


(6.27 + 0.14) x 10~ n 


(6.49 + 0.15) x 10-" 



a) The photon fluxes Fo_i and Fo,2 are expressed in phcm -2 s _1 keV~'at £0 = 1 keV. 

b) Error estimates (lcr) in the fit parameters are obtained by fixing all parameters except the parameter 
of interest at their optimum values and subsequently determining, where the parameter of interest 
reaches a A^?_ t step of 7.038 from the global minimum. This increment corresponds to a lcr step 
for 6- 1 degrees of freedom. 



We used logparabolic functions in order for the spectra 
to bend downwards to the optical regime, as required by the 
broad-band spectrum (se e e.g. the s pectrum from the first 
broad-band campaign by Iden Hartog et al.1 120071) . If a tradi- 
tional model, a combination of a black body and a power law, 
was chosen, the soft power-law component would keep in- 
creasing towards lower energies, which can only be compen- 
sated by increasing the absorption column. T he latter would be 
forced to unrealistically high values (see e.g. iRea et alj|2007al 



who derived N H = (1.00 ± 0.01) x 10 zz cm' 2 ). Fixing the 



A'h to the value published by iDurant & van Kerkwiikl (2006b, 
A^h = (0.64 ± 0.07) x 10 22 cirT 2 ) from an accurate analysis 
of data of the Reflection Grating Spectrometer aboard XMM- 
Newton yielded, however, also poor fits at low energies be- 
tween ~0.6 keV and ~1 keV in both observations B and C. 
Therefore, we treated Nh as a free parameter. For both observa- 
tions we find an absorption column of (0.57 ±0.02) x 1 22 cm -2 , 
which is so mewhat lower than, b ut consistent with the value ob- 
tained by Durant & van Kerk wiik (2006b). Therefore, assum- 
ing that the origin of the absorption is Galactic and thus con- 
stant, we fix JV H to 0.57 x 10 22 cm" 2 fo r the rest of the XMM- 
Newton analyses. We note that lGiiver et al. (2008) find similar 
values for Nh by applying a surface thermal-emission and mag- 
netosphere scattering model to the same XMM-Newton data in 
a parallel analysis. 

Next, we consider the total XMM-Newton spectral energy 
window up to ~12 keV. An additional model component is re- 
quired to describe at least a part (< 20%) of the harder (~5- 
12 keV) X-ray emission from both observations B and C. For 
this component a second logparabolic function is introduced, 
instead of the traditionally introduced power-law component. 
Furthermore, a third logparabolic component with fixed model 
parameters prescribed by the ISGRI/SPI joint fit above 20 keV 
(see Sect. 13. 1.2b contributes in this energy band and is included 
in the total fit. Without the need to introduce an XMM-Newton- 
ISGRI inter-calibration factor, a good fit was obtained to the 



high-energy parts of the XMM-Newton spectra. In Table|5]the 
fit results are presented for the three XMM-Newton observa- 
tions. 



However, for both long observations (obs. B and C) sta- 
tistically unacceptable fits are obtained employing this empir- 
ical approach. Namely, significant residuals are apparent be- 
low 2 keV, because the logparabolae cannot follow the spectral 
shape exactly at these low energies. The small deviations from 
the optimum fit (hardly visible in Fig. [3]) are of the order of a 
few percent for energies below 2 keV. We use the spectral fits, 
however, primarily for the conversion of absorbed fluxes into 
unabsorbed fluxes (data points shown in Fig. [3} by calculating 
for each energy bin the fraction of energy which is absorbed. In 
Sect. l3.3l we use these high-statistics unabsorbed total flux mea- 
surements in even broader bins for deriving the pulsed fraction 
as a function of energy, and the uncertainty in that analysis is 
dictated by the errors in the fluxes from the phase resolved anal- 
ysis. Therefore, the empirical two plus one (fixed) logparabolae 
fits assuming a constant absorption column of 0.57 x 10 22 cm" 2 
are sufficiently accurate for our purposes. 



Note from Table|5]that the total unabsorbed 2-10 keV en- 
ergy flux is consistent with being constant during the time span 
of th e three XMM-New ton observations. Th i s is in agreement 



with Rea et al. 



Gonzalez et al. 



1 2007ah and iGonzalez etail J2008h . however, 
d2008l) report a flux increase of (15±3)% (2- 
10 keV) between the epochs of our XMM-Newton observa- 
tions A, B and C and the epochs after the bursting activities 
started in 2006. The flux levels around 10 keV of the three 
XMM-Newton spectra are rather similar. Fig. [3] showing only 
the XMM-Newton obs. B spectrum, is therefore a good repre- 
sentation of the total high-energy (>0.5 keV) spectrum, partic- 
ularly showing the abrupt transition from soft to hard X-rays. 
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Fig. 4. 20-160 keV ISGRI pulse profile of 4U 0142+61. The 
significance is 6.5<x using a Z^ test (fit shown as solid curve). 
The fitted DC level is indicated with a horizontal line. Also 
indicated are the phase intervals for Ph I, II and III (see Tabled 
with vertical grey lines and with colours. 



3.2. Timing analysis: pulse profiles 

In this section we first present the time-averaged INTEGRAL 
IBIS-ISGRI pulse profiles of 4U 0142+61. Then, we will ex- 
amine the XMM-Newton EPIC-PN and ASCA-GIS pulse pro- 
files in single observations and test these for variability with 
time and energy, respectively. Finally, we compare the time- 
averaged ISGRI pulse profiles with the time-averaged XMM- 
Newton pulse profiles. 



3.2.1. INTEGRAL ISGRI pulse profiles 



Using the procedure described in Sect. I2.1.2I it is possible to 
detect significant pulsed emission from 4U 0142+61 in the 
ISGRI data. For the energy range 20-160 keV a maximum sig- 
nificance is reached of 6.5<x using a X\ test. The 20-160 keV 
pulse profile, shown in Fig.H] is broad and appears to be dou- 
ble peaked with a DC level in a narrow phase interval of width 
-0.2. 

Pulse profiles for the differential energy intervals 20- 
50 keV and 50-160 keV are shown in Fig. |7]E, F. The signifi- 
cances of the ISGRI pulse profiles are A. Act and 4.7<x, respec- 
tively, and are both more significant than the RXTE-HEXTE 
profiles f or similar energy ran ges (3.4cr and 2. Oct, respectively) 
shown in iKuiper et alj (120061) . The profile for the 20-50 keV 
band (Fig. 0E) shows the two peaks separated ~0.5 in phase. 
The 50-160 keV profile (Fig. 0F), however, shows only a 
shoulder to the first pulse at the location of the second pulse 
(near phase ~0.5) in the 20-50 keV profile. The minimum (DC 
level) is in both cases located in the same narrow phase interval. 



In order to quantify a possible pulse-morphology change 
as a function of energy, we performed a statistical method 
based on a combination of a Pearson^ 2 (see e.g. Sect. 11.2 
of Eadie et al. 1971) and a Run test (see e.g. Sect. 11.3.1 of 
Eadie et alJll97ll) . comparing the shapes of the (independent) 
20-50 keV and 50-160 keV pulse profiles. We did not apply 
the Kolmogorov-Smirnov test, because these background dom- 
inated (binned) ISGRI pulse phase distributions have different 
backgrounds and DC levels. In our approach we first obtained 
an appropriate model function describing the shape of the 20- 
50 keV pulse profile by fitting a truncated Fourier series using 
three harmonics. Next, we fitted the 50-160 keV distribution 
(in 20 bins) in terms of a constant and a function with free 
scale describing the shape of the 20-50 keV pulse profile (dof 
= 18). We investigated not only the absolute values of the de- 
viation of the data from the optimum fit, but also the signs of 
the deviations by combining the probabilities from the inde- 
pendent Pearson^- 2 test (^V) and Run test (P rm ) in an overall 
joint probabili ty, P = P r i ■ P nm (l - lnP^ • P run ) (see e.g. Sect. 
11.6 of lEadie et al.lll97ll) . In our case the joint probability is 



0.104, thus we cannot claim a significant difference in shape 
between the 20-50 keV and 50-160 keV ISGRI pulse profiles. 

3.2.2. XMM-Newton and ASCA pulse profiles: time 
variability 

The ASCA-GIS an d RXTE-PCA pulse profiles published by 
Kuiper et alj ((2006) suggested that the pulse profiles below 



~2-3 keV are significantly different. However, the RXTE-PCA 
pulse profiles were averages over years and the RXT E-PCA is 



not sensitive below 2 keV, contrary to ASCA-GIS . Dib et al. 



(2007) used the 10 years of RXTE monitoring above 2 keV to 
show long-term variability in the pulse-profile shape. However, 
the reported effect is small; mainly the dip of emission be- 
tween the two peaks got shallower between 2002 and 2006. 
Therefore, we used the ASCA-GIS and XMM-Newton obser- 
vations to create pulse profiles in the same differential energy 
bands to test for time variability. The differences in the en- 
ergy resolutions and responses between the two instruments 
has negligible impact on the comparison of the profile shapes 
in the selected broad energy intervals. 

Fig.|5]shows ASCA-GIS and XMM-Newton pulse profiles 
in the energy range 0.8-2.0 keV. This figure strongly suggests 
that the pulse profile taken with ASCA-GIS is different from 
the three XMM-Newton pulse profiles. For ASCA-GIS, the 
pulse that peaks around phase 0.6 is higher than the pulse that 
peaks around phase 1.1, while it is the other way around for 
all three XMM-Newton profiles. The latter profiles seem rather 
similar. We tested this applying the combination of the Pearson 
X 1 and Run tests described in Sect. 13.2.11 The XMM-Newton 
profile with the highest statistics (Obs. B) was used as tem- 
plate. XMM-Newton profile A, with the lowest statistics, differs 
from profile B at the 1.8cr level, and profile C differs from pro- 
file B at the 2.6<x level. Slight differences in shapes are visible 
around the peaks of the two pulses, but we do not consider this 
evidence for significant time variability between the XMM- 
Newton observations. To the contrary, the ASCA-GIS profile 
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Fig. 5. 4U 0142+61 pulse profiles for 0.8-2.0 keV of ASCA- Fig. 6. 4U 0142+61 pulse profiles for 2.0-8.0 keV of ASCA- 
GIS and the three XMM-Newton observations (A, B, C, GIS and the three XMM-Newton observations (A, B, C, 
Table [3}. Phase intervals indicated as in Fig. [4] The profile TableQ. The relative normalization scales approximately with 
of XMM-Newton observation C required a shift in phase by the total pulsed flux (2-8 keV; see e.g. Table 8). See also the 
-0.028 for alignment with the profile of observation B (con- caption of Fig. [5] 
sistent with the other profiles). The relative normalization is 
consistent with the total pulsed flux (0.8-2.0 keV) being ap- 
proximately constant (see Table|8]l. 



for 0.8-2.0 keV differs from XMM-Newton profile B at the 
16.6<x level. The above-mentioned differences in peak heights 
in the ASCA-GIS and XMM-Newton profiles become also 
clear when comparing the fluxes in these pulses (see Sect. 13.3.21 
and Table |H). The pulse peaking at phase 0.6 contains ~18% 
more flux than the average pulse during the XMM-Newton ob- 
servations. The pulse peaking at phase 1.1 contains ~34% less 
flux than the average pulse during the XMM-Newton observa- 
tions. 

Similar conclusions can be reached for the energy band 
above 2 keV for which the profiles are shown in Fig. [6] XMM- 
Newton profiles A and C differ from profile B at the 1.3cr 
and 3.5cr levels, respectivehjj; the ASCA-GIS profile for 2.0- 
8.0 keV deviates at the 9.3<x level. The variation is more drastic 
than the difference in profile shape may indicate. The flux in the 
GIS pulse peaking at phase 1 . 1 is 44% lower than the averaged 
flux in the same pulse in the XMM-Newton profiles. The flux 
in the GIS pulse peaking at phase 0.6 is only 12% reduced in 
flux. The difference is evident in Fig. [6] in which the profiles 
are approximately scaled to the measured 2-8 keV total pulsed 
flux. 

We conclude that we found significant variability in the 
overall pulse-profile shapes for energies below 8 keV from 



3 iRea et al.l ( 12007 ah . analysing the same XMM-Newton observa- 
tions obtained in 2004, also concluded that there are no significant 
differences in the profile shapes. 



the epoch of the ASCA observations (July /August 1999) to 
the epochs of the XMM-Newton observations in the years 
2003-2004 (see Table |3), Particularly in the relative fluxes 
and spectra of the two pulses the change in shape is more 
dra stic than r eporte d before from RXTE-PCA observations 
by iDib et all d2007l) and from more XMM-Newton observa- 
tions by Gon zalez et al.l d2008h . iMorii et all d2005bl) also re- 
port pulse-shape changes from the ASCA-GIS data. They use 
the 1999 observations (which we also use here) as a template 
to test two shorter observations from 1994 and 1998. Due to 
the low statistical quality of these data sets only significant 
changes below 3 keV could be claimed. IDib et al. ( 2007b no 



tice that the 1998 ASCA-GIS profile presented by IMorii et al. 



d2005bl) agrees with the pre-gap RXTE profile while the 1999 
one does not. Interesti ngly, the profile measured with Chandra 
dGonzalez et al.ll2008l) in May 2000 seems to be a transition 
between the ASCA-GIS (1999) and XMM-Newton profiles. If 
the deviant ASCA-GIS profile is c aused by the glitch prior 
to this observation as proposed by IMorii et al.l d2005al) . then 
it took about a year to get to the stable situation as measured 
with XMM-Newton i n 2003-2004 an d as can be seen from the 
RXTE pulse profiles dDib et al.ll2007b . 



Apparently, some relatively fast reconfiguration in the 
pulse shape is ongoing during the 1999 ASCA-GIS and 2000 
Chandra observatio ns, which might be due to a proposed glitch 
dMorii et al.ll2005bl) or other drastic event. 
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Fig. 7. 4U 0142+61 pulse profiles from soft to hard X-rays. 
In panels A, B and C, pulse profiles in the energy ranges 0.8- 
2.0 keV, 2.0-4.0 keV and 4.0-8.0 keV are shown which are 
the sums of the XMM-Newton observations A, B, and C. Panel 
D presents the RXTE-P CA pulse profile in the energy band 
8.0-16.3 keV taken from lKuiper et alj (1200 6). Finally, in pan- 
els E and F INTEGRAL pulse profiles are shown in the energy 
ranges 20-50 keV and 50-160 keV. Phase intervals are indi- 
cated as in Fig. [4] 



component that is apparent in all energy bands is the pulse that 
peaks around phase 1.1. Whereas the pulse that peaks at phase 
0.7 in the XMM-Newton band is vanishing and no longer vis- 
ible in the RXTE-PCA band, in the INTEGRAL band a pulse 
is visible around phase 0.5. In Sect. [33] we will disentangle the 
contributions from different components using phase-resolved 
spectroscopy. 

3.3. Timing analysis: total-pulsed and phase-resolved 
pulsed spectra 

Total-pulsed spectra were created by the procedure described in 
Sect. 12.1.21 An example of how the excess counts in the pulsed 
components were extracted above a flat background level is 
shown in Fig. [4] 

We performed phase-resolve d spectroscopy fo r the first 
time for just the pulsed emission. iRea et al.l d2007al) also per- 
formed phase-resolved spectroscopy with the XMM-Newton 
data. However, they did not separate the pulsed component 
from the dominant (on average ~85%) DC component, which 
has a different spectrum and is invariant with phase. For our 
phase-resolved spectroscopy, we selected three phase intervals 
for which we want to extract high-energy spectra, using the 
morphology of the ISGRI 20-160 keV pulse profile (Fig. gj 
for defining the phase boundaries: Ph I contains the smaller 2nd 
pulse in the INTEGRAL band, Phil the DC level and Phlll 
contains the main pulse at hard X-rays, as well as the main 
pulse for energies below 10 keV. The definitions of the phase 
intervals are given in Tableland are indicated with grey verti- 
cal lines in all figures with pulse profiles. 

Table 6. Selected phase intervals for extraction of high-energy 
spectra, using the pulse-shape morphology of the INTEGRAL- 
IBIS-ISGRI 20-160 keV pulse profile (Fig.©. 



Phase interval 



Component 



Phi [0.35,0.666) 

Phil [0.666,0.851) 

Phlll [0.0,0.35) V [0.851,1.0) 



Secondary INTEGRAL pulse 
DC level INTEGRAL profile 
Main INTEGRAL pulse 



3.2.3. Pulse-profile changes with energy 

Kuiper et al] (12006) presented pulse profiles as a function of 
energy from ASCA-GIS, RXTE-PCA and HEXTE . From that 
work an d e.g. Ilsrael et all dl999h : bib et al.l d2007l) : IRea et al.l 
(2007a) it is clear that the pulse profiles significantly change 
morphology with energy. In Sect. I3.2.2l we have shown that the 
pulse profiles of the 2003-2004 XMM-Newton observations 
are statistically the same. Therefore, we use in this section the 
summed 'time-averaged' pulse profile using all three XMM- 
Newton observations for comparison with maximal statistics 
with the profiles derived at higher energies with RXTE-PCA, 
and INTEGRAL. 

In Fig. Qtime-averaged pulse profiles from XMM-Newton, 
RXTE-PCA, and ISGRI, are shown in 6 differential energy 
bands between 0.8 keV and 160 keV. It is evident that the mor- 
phology changes are ongoing throughout the energy range. One 



3.3.1. INTEGRAL and RXTE 

For the INTEGRAL data we constructed four broad-band pulse 
profiles in the energy ranges 20^-0 keV; 40-80 keV; 80- 
160 keV and 160-300 keV. From these pulse profiles we got 
three total-pulsed flux measurements and one upper limit. 

For RXTE-PCA, we extracted 14 pulse profiles in the en- 
ergy band 2.5-31.5 keV for which we determined the number 
of excess counts (see Sect. l2.21 >. The last energy bin (23.8-31.5 
keV) yielded no significant pulse profile, but the INTEGRAL 
and RXTE-PCA spectra nicely bridge the minimum luminosity. 

The time-averaged total-pulsed spectra derived from 
INTEGRAL and RXTE-PCA data are presented in Fig. g] The 
spectra are fitted simultaneously using a logparabola for the 
soft X-rays below 10 keV and either a power law or a log- 
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Fig. 8. In the top panel the phase-resolved pulsed-emission 
spectra and fits of INTEGRAL and RXTE-PCA are presented. 
The data points with a square symbol are INTEGRAL mea- 
surements and those without a marker are from RXTE-PCA. 
In black is presented the total-pulsed spectrum; in blue Phi; 
in green Ph II; and in magenta Ph III. The arrows indicate the 
flux measurements with a significance less than 1 .5<x. For com- 
parison the INTEGRAL/XMM-/Vewfo«-B total-spectrum fit is 
drawn in black. Finally, in grey is plotted the power-law fit for 
the INTEGRAL total spectrum and the power-law fit for the 
total-pulsed spectrum. In the bottom panel 77 is presented in the 
same colour scheme as the spectra. 77 is defined as the fraction 
of the pulsed emission in a selected phase interval Ph I, Ph II or 
Phlll of the total-pulsed emission, i.e. the sum equals unity. 

parabola for the hard X-rays. The fit results are summarized 
in Table [7] In Fig. [8] the fit with two logparabolae is plotted 
in black while the logparabola and a power-law fit is drawn 
in grey. Both fits to the INTEGRAL data are fully acceptable, 
but the single power-law model uses less parameters. It can be 
seen that if the hard X-ray spectrum has a power-law shape 



Table 7. RXTE-PCA and INTEGRAL spectral-fit parameters 
for the total-pulsed (TP) spectrum and for three phase inter- 
vals (see Table©. The N H is fixed to 0.57 x 10 22 cm" 2 . The 
subscripts below the fit parameters indicate if this parameter 
accounts for the soft (s, <10 keV) or hard (h, >10 keV) X- 
ray part of the spectrum. The normalizations (Fq s and Foh) are 
taken at Eq - 1 keV and have the units phcirT 2 s keV -1 . 



TP 


Fit 1 


Fit 2 


#s 


2.27 ±1.00 


2.82 ±1.08 


A 


1.71 ±0.90 


1.10 ±0.96 


Fo. 


(1.34 ± 0.38) x 10- 2 


(1.79 ± 1.23) x lO" 2 


r h 


0.40 ±0.15 




<>'h 




-2.83 ±2.16 


A 




1.00 ±0.67 


-foh 


(9.17 ± 4.47) x 10- 6 


(0.34 ± 2.62) x lO" 8 


xl (dof) 


0.78 (13) 


0.50(12) 


ft-iokev (erg cm -2 s" 


») (6.82 ± 1.77) x lO" 12 


(6.85 ± 0.93) x 10-' 2 


-F20-150keV 


(2.68 ± 1.34) x 10-" 


(2.77 ± 0.68) x 10-" 


Phi 


a s 


2.90 ±2.14 


2.89 ±3.66 


A 


2.56 ±2.88 


2.42 ±3.35 


F 0s 


(0.67 ± 1.71) x lO" 2 


(0.63 ± 1.74) x 10- 2 


r„ 


0.12 ±0.20 




<>'h 




-4.00 ± 4.62 


A 




1.31 ±0.72 


Foh 


(5.78 ± 8.81) x 10~ 7 


(0.67 ± 5.47) x 10- 9 


X' (dof) 


1.43(13) 


0.69 (12) 


^2-iokeV (ergcirT 2 s~ 


') (1.62 ± 0.60) x lO" 12 


(1.14 ± 3.64) x 10~ 12 


-f20-150keV 


(5.99 ± 3.29) x 10- 12 


(6.89 ± 6.42) x 10-" 


Phil 


<!"s 


2.95 ± 1.71 




/3s 


1.59 ± 1.56 




F 0s 


(5.25 ± 0.72) x lO" 3 




Xl (dof) 


0.86 (15) 




^2-iokeV (erg cm" 2 s" 


') (1.30 ± 0.03) x 10~ 12 




Phlll 


<!"s 


1.63 ±0.59 


2.68 ± 0.79 


/3s 


2.07 ± 0.44 


1.01 ±0.70 


F 0s 


(5.55 ± 2.34) x 10- 3 


(1.02 ± 0.52) x lO" 2 


r„ 


0.50 ±0.11 




<>'h 




-2.99 ± 1.73 


A 




1.09 ±0.54 


-F()h 


(1.11 ± 0.39) x 10~ 5 


(0.26 ± 1.67) x 10~ 7 


Xl (dof) 


1.1(13) 


0.65 (12) 


ft-iokev (ergcirrV 


') (4.72 ± 0.52) x lO" 12 


(4.82 ± 0.47) x lO" 12 


-f20-150keV 


(2.06 ± 0.82) x 10 -11 


(2.07 ± 1.09) x 10-" 



with r = 0.4 that the pulsed fraction could become as high as 
100% around ~250 keV. If the spectral shape is logparabolic, 
the pulsed fraction would not increase and it would stabilise 
around ~30%-4Q% 

For the three differential phase intervals (see Table |6]l also 
spectra are extracted and plotted in Fig. [8] Immediately clear is 
that Ph III contains the bulk of the energy of the pulsed emis- 
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Fig. 9. For 4 (soft) X-ray and INTEGRAL observations we present the total, total-pulsed and phase-resolved pulsed spectra (top 
panels) and 77 (bottom panels, defined as the fraction of the total pulsed emission). The (soft) X-ray data from left to right are 
ASCA GIS, XMM-Newton PN obs. A, B and C. All data are equally binned for easy comparisons. For the total spectra, the best 
fit is shown. In the top panel for ASCA-GIS the best fit to the total spectrum of XMM-Newton obs. B and INTEGRAL are shown 
for clarity. The total spectra and total-pulsed spectra are plotted in black. The phase resolved spectra for Ph I, Ph II and Ph III (see 
Tabled are plotted in blue, green and magenta, respectively. The INTEGRAL phase-resolved spectra are indicated in each figure 
with a square symbol. The data points with arrows indicate that the flux values have significances less than 1.5cr. The positive 
lcr error is drawn. In the bottom panels 77 is presented in the same colour scheme as the spectra. 



sion. The spectrum shows a similar transition (from soft to 
hard) as the total-pulsed spectrum. Also visible is that the soft 
X-ray spectrum (<10 keV) of Phlll is quite a bit harder than 
the soft X-ray spectra of Phi and Phil. While below 10 keV 
Ph I and Ph II form one pulse, in the INTEGRAL band only in 
Phi is a pulse visible. From the spectrum and fit in Fig.[8]it can 
be seen that this transition for Ph I is extremely drastic from the 
very soft to the very hard. 

We can not present an accurate pulsed fraction using 
RXTE-PCA data. The pulsed fraction is defined as the pulsed 
flux (determined as described in Sect. |2j divided by the to- 
tal flux from the point source (pulsed + DC). The RXTE- 
PCA is a non-imaging instrument, making it difficult to de- 
termine a reliable total flux and spectrum due to confusion 
with RX J0146.9+ 6121, a closeby high-mass X-ray binary 
(Motchet al. 1991), and underlying Galactic -ridge emission 
(Valinia & Marshall 1998). Alternatively, we present 77 in the 
bottom panel of Fig. [8] 77 is defined as the fraction of the pulsed 
emission in a selected phase interval of the total-pulsed emis- 
sion. From the 77 of Ph III (magenta) it can be seen that the 
contribution of this phase interval to the total-pulsed emission 
increases from ~65% around 3 keV to practically 100% around 
10 keV. The contributions of Ph I and Ph II gradually decrease 
from ~20% at 2 keV to (consistent with) no contribution at 
10 keV. The contribution of Phil remains consistent with 0% 
in the INTEGRAL band (DC level), but the contribution to the 



pulsed emission from Phi increases again as a result of the 
small 2nd peak in the INTEGRAL band. Because the pulsed 
contributions in both phase intervals Phi and Phil similarly 
decrease to ~Q% at 10 keV, and the pulse shape in the soft X- 
ray band seems to be of one component, we infer that the small 
2nd pulse in the INTEGRAL band is unrelated to the soft X- 
ray pulse. It likely originates from a different site (height) in 
the magnetosphere but appears around the same phase as the 
soft X-ray pulse. 

It is important to remember that the time-averaged 
RXTE-PCA spectra are assembled between March 1996 and 
September 2003, while the INTEGRAL observations are taken 
from December 2003 to August 2007. It is then remarkable 
how well the spectra (total pulsed and Ph III) connect in Fig. [8] 
suggesting stable emission (geometry) above ~10 keV over 1 1 
years. 

3.3.2. XMM-Newton and ASCA 

For the three XMM-Newton and ASCA observations we also 
determined the total-pulsed spectra as well as the pulsed emis- 
sion spectra for the three selected phase intervals. Systematic 
uncertainties in the response for the timing mode of XMM- 
Newton above ~8 keV prevented us from using the higher- 
energy data. For this reason the observational gap between 
XMM-Newton and INTEGRAL is about 12 keV. 
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Table 8. ASCA-GIS and XMM-Newton (A, B & C) fits for the total-pulsed (TP) spectra and for spectra of three phase intervals 
(Phi, II and III, see Tabled for4U 0142+61. The A^h is fixed to 0.57 x 10 22 cm 2 . For each observation the three parameters of 
the parabola are given. The normalizations (Fq) are taken at Eq — 1 keV and have the units phcirr 2 s" 1 keV" 1 . 



TP 




ASCA-GIS 


A 


B 


C 


(i- 




1.88 + 0.30 


2.19 + 0.38 


1.24 + 0.15 


1.12 + 0.17 


P 




2.03 ± 0.48 


0.65 ± 0.55 


2.43 ± 0.26 


2.41 ± 0.27 


F 




(8.49 ± 0.77) x 10" 3 


(7.93 + 0.10) x 10" 3 


(8.17 + 0.41) x 10" 3 


(7.98 ± 0.45) x 10" 3 


x\ (dof) 




1.8 (7) 


0.29 (6) 


0.88(15) 


0.86(15) 


^o.8-2 kcv (ergcm"V 


') 


(11.47 + 0.70) x 10" 12 


(10.96+ 1.38) x 10" 12 


(12.64 + 0.50) x 10" 12 


(12.72 + 0.56) x 10" 12 


^2-iokeV (erg cm" 2 s" 


l ) 


(5.55 ± 0.33) x 10" 12 


(9.02+ 1.20) xlO" 12 


(8.16 + 0.35) xlO" 12 


(9.17 + 0.39) xlO" 12 


Phl 


a 




1.54 + 0.27 


2.72 ± 0.47 


1.57 + 0.20 


1.59 + 0.20 


fi 




3.76 + 0.68 


0.39 ± 0.84 


2.71+0.45 


2.79 ± 0.45 


F 




(3.49 ± 0.28) x 10" 3 


(2.57 ± 0.38) x 10" 3 


(2.53 + 0.14) x 10" 3 


(2.57 + 0.14) x 10" 3 


x\ (dof) 




0.66 (7) 


0.30 (6) 


0.61 (15) 


0.50(15) 


^o.8-2 kuv (erg cm" 2 s" 


l ) 


(4.58 ± 0.46) x 10" 12 


(3.29 ± 0.80) x 10" 12 


(3.56 + 0.27) x 10" 12 


(3.37 + 0.25) x 10" 12 


^2-iokeV (erg cm" 2 s" 


') 


(1.39 ± 0.14) x 10 -12 


(1.83 + 0.52) x 10" 12 


(1.50 ±0.16) xlO" 12 


(1.36 + 0.12) x 10~ 12 


Phil 


a 




1.51+0.31 


2.24 ± 0.39 


0.94 + 0.18 


0.97 + 0.17 


P 




2.81+0.52 


0.83 ± 0.63 


3.27 ± 0.34 


2.97 + 0.31 


F 




(1.79 + 0.16) x 10" 3 


(1.80+ 0.23) x 10" 3 


(1.75 + 0.10) x 10" 3 


(1.51 + 0.08) x 10" 3 


Xl (dof) 




1.73 (7) 


1.09(6) 


0.96(15) 


1.35(15) 


Fo.8-2keV (erg cm" 2 s" 


') 


(2.50 + 0.26) x 10" 12 


(2.44 + 0.50) x 10" 12 


(2.77 ± 0.20) x 10" 12 


(2.41+ 0.18) x 10" 12 


^2-iokcV (erg cm" 2 s" 


') 


(1.13 + 0.10) xlO" 12 


(1.65 + 0.39) xlO" 12 


(1.51 ±0.10) xlO" 12 


(1.38 + 0.10) x 10" 12 


Phlll 


a 




1.87 + 0.33 


1.74 + 0.43 


1.05 + 0.14 


0.95 + 0.14 


P 




1.36 + 0.44 


0.95 ± 0.55 


2.31+0.21 


2.41 ± 0.20 


F 




(3.03 ± 0.37) x 10" 3 


(3.42 ± 0.55) x 10" 3 


(3.86 + 0.19) x 10" 3 


(3.98 + 0.19) x 10" 3 


Xl (dof) 




1.72 (7) 


0.27 (6) 


1.01 (15) 


1.54(15) 


^o.8-2 kcv (erg cm" 2 s" 


') 


(4.23 + 0.58) x 10" 12 


(5.09+ 1.17) x 10" 12 


(6.30 + 0.37) x 10" 12 


(6.65 + 0.36) x 10" 12 


^2-iokcV (erg cm" 2 s" 


') 


(3.01 + 0.28) x 10" 12 


(5.37 ± 0.82) x 10" 12 


(5.13 + 0.24) x 10" 12 


(5.63 + 0.22) x 10" 12 



In the top panels of Fig. [9] the results of ASCA-GIS and 
XMM-Newton obs. A, B and C are presented together with the 
INTEGRAL results of Fig. [8] The colour scheme is the same 
as in Fig. [8] In black are drawn the total-pulsed spectra and 
(except for ASCA-GIS) the (binned) total spectra. To the to- 
tal spectra of XMM-Newton and INTEGRAL the best fits are 
drawn. In the panel for ASCA-GIS the fit to the total spectrum 
for XMM-Newton obs. C and INTEGRAL is drawn. 

We have fitted each spectrum with a single logparabolic 
function. These fit results are given in Table [8] The parame- 
ters of a logparabolic function are not easy to compare, but the 
flux values indicate that there are some differences between the 
different observations, consistent with our findings in the com- 
parison of the pulse shapes. 

Most remarkable are flux differences between the ASCA- 
GIS observation and the XMM-Newto n ones. (The latter ar e 



consistent with the fluxes reported by iGonzalez et al. 2008) 



While the 0.8-2 keV total-pulsed fluxes of all observations are 
comparable, the 2-10 keV total-pulsed flux of ASCA-GIS is 
significantly lower than the rest. From the phase-resolved spec- 
tra it becomes clear that this difference is due to the signif- 
icantly lower flux (2-10 keV) in the main pulse (Phlll). In 



this phase interval the ASCA-GIS 2-10 keV flux is less than 
in the XMM-Newton observations, while the 2-10 keV flux for 
phases Phi and Phil are slightly lower than, or comparable 
to those measured with XMM-Newton. We note that the 0.8- 
2 keV ASCA-GIS flux in Ph III is somewhat lower than in the 
other observations. Note also that only for ASCA-GIS the 0.8- 
2 keV flux in Ph I is higher than, or comparable to, this flux in 
Phlll. 

Table [8] shows that over the three XMM-Newton observa- 
tions the total-pulsed flux in the 0.8-2 keV and 2-8 keV bands 
as well as the fluxes in all three phase intervals are statisti- 
cally stable. The latter is consistent with the conclusion that 
the shape of the pulse profiles did not change measurably in 
this time interval. 

3.3.3. Comparison 

To determine the pulsed fraction as a function of energy we 
have binned the total spectra in the same binning as the total- 
pulsed and phase-resolved pulsed spectra. The pulsed fraction 
is defined as the pulsed flux (determined as described in Sect.|2]i 
divided by the total flux from the point source (pulsed + DC). 
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For the instruments ASCA-GIS and RXTE-PCA we do not 
have accurate total spectra, therefore we also consider 77. 

We start with the pulsed fraction. In the previous sections 
we showed that XMM-Newton observations A, B and C exhibit 
statistically identical pulse shapes and pulsed fluxes. Therefore, 
we have combined these observations by taking weighted av- 
erages. In Fig. [TO] the pulsed fraction as a function of energy 
is presented for the averaged XMM-Newton and INTEGRAL 
observations (black data points). Also presented in this fig- 
ure are the pulsed fractions for the three selected phase inter- 
vals (coloured data points). It can be seen in this figure that 
the pulsed fraction in the XMM-Newton band increases from 
~10% around 1 keV to ~20% around 8 keV. This trend con- 
tinues in the INTEGRAL band where the pulsed fraction is 
increasing to p ossibly as high as ~40% in the 40-80 keV band. 
Rea et al.1 d2007al) presented the pulsed fraction of the funda- 
mental and first harmonic sine functions for the XMM-Newton 
pulse profiles as a function of energy below 10 keV. These are 
consistent with the total pulsed fractions we show in Fig.fTOl 

Fig. [10] also shows that the trend set by the total-pulsed 
fraction is followed closely by the pulsed fraction that can be 
assigned to Ph III (magenta), as we have already discussed in 
Sect. [33J] Phlll alone accounts for ~5%-12% of the pulsed 
fraction between 1-8 keV. The other two phases do not show 
the same trend. Ph I, the 2nd pulse in the INTEGRAL band, ac- 
counts for ~5% at 1 keV decreasing to ~3% at 8 keV. Finally, 
Ph II is about constant in this energy range at a level of ~2.5%. 
These trends are also visible in the phase-resolved pulsed spec- 
tra in Fig. [9] showing that at higher energies around 8 keV 
Ph III is dominating over the other two phase intervals. 

To investigate this further, the distribution of 77 is even more 
illustrative, as can be seen in Fig.[9]and Fig.QT] the contribution 
of Ph III to the total pulsed emission is steadily increasing with 
energy, reaching a fraction of ~65% at 8 keV. For the avarage of 
the above reported RXTE-PCA observations (Fig.[8]the contri- 
bution from Phlll increases even to ~100% around 10 keV. In 
the INTEGRAL band this contribution is again down to ~60%- 
70%, due to the new pulse in Ph I. 

In Fig. [9] is also clearly visible that during the ASCA-GIS 
observation the 77 of Phlll below 2 keV is less than that of 
Phi. In all observations the 77 of Phi is decreasing with en- 
ergy from ~30%-40% down to <5% around 8 keV. However, 
above 10 keV the contribution of Ph I is increasing again due to 
the 2nd pulse in the INTEGRAL band. Ph II, which was chosen 
as the DC level of unpulsed emission in the INTEGRAL band, 
shows first a constant 77 of ~20% and it starts decreasing from 
2 keV gradually to 0% in the INTEGRAL band. 

4. Summary 

In this paper we have presented new and more detailed char- 
acteristics of 4U 0142+61 in the hard X-ray regime (>10 keV) 
and studied these in a broader high-energy view using archival 
soft X-ray data (<10 keV) in consistent analyses. Many new 
or more accurate results are reported which should be consid- 
ered in theoretical modelling of the high-energy emission from 
AXPs. Particularly new for studies of AXPs is the performed 
phase resolved spectroscopy over the total high-energy band, 




Energy [keV] 

Fig. 10. Pulsed fraction as function of energy. Below 10 keV 
the averaged XMM-Newton pulsed fraction of observations A, 
B and C are plotted. Above 10 keV the INTEGRAL pulsed 
fraction is plotted. Consistent colour coding as Fig[8]and Fig. [9] 
Total-pulsed emission (black), Ph I (blue), Ph II (green), Ph III 
(magenta). 



revealing distinctly different components contributing to the to- 
tal emission. We will first present a summary of the results. 

4.1. Total high-energy emission of 4U 0142+61 

1) ISGRI measures with an effective exposure of 2.4 Ms the 
most accurate time-averaged total spectrum between 20 and 
230 keV which is fully consistent with a power-law shape with 
index 0.93+0.06 (Sect. UlU Fig.\3\and Tabled- 

2) Combining ISGRI and SPI spectral data provides the first 
evidence at the ~4<x level for a spectral break with maximum 
luminosity at an energy of ~270 keV, assuming a logparabolic 
shape of the hard X-ray spectrum (Sect. \3. L2\ and Fis. lit. 

3) The total high-energy spectrum, combining spectra mea- 
sured with XMM-./Vewfo77 and ISGRI, can be described sat- 
isfactorily with the sum of three logparabolic functions, two 
dominating the spectru m for ene rgies below 8 keV and one 
above 20 keV (see also lRea et al.l l2007b). Minimum values for 
the luminosity are reached around 10 keV {Sect. \3.1.4\ Fig. \3\ 
and Table\5^ . 



4.2. Long-term variability in the total emission of 
4U 0142+61 

1) ISGRI does not see evidence for significant long-term vari- 
ability of the source flux and power-law spectral index on one- 
year time scales between December 2003 and August 2006 at 
the 17% (lo-) level (Sect. UJl\ Fig.U\&®and Table^. 
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Fig. 11. Contribution to the total-pulsed emission rj of the three 
selected phase intervals from the averaged XMM-Newton rj 
of observations A, B and C and for INTEGRAL (>10 keV). 
Colour coding consistent with Fig.fTOl 



2) Two short INTEGRAL ToO observations following burst- 
ing behaviours of 4U 0142+61 measured with RXTE show the 
AXP to be in a state consistent with the time averaged flux 
and spectral shape. The total flux was stable within ~20% (lcr; 
Sect. UlH . 

4.3. Pulse profiles 

1) ISGRI measures above 20 keV significant hard-X-ray pulse 
profiles up to 160keV (Sect. \3.2J\ and Fis. l4ii. 

2) The morphologies of the profiles measured below 10 keV 
with XMM-JVewto« and ASCA-GIS differ from the hard-X-ray 
profile; one of the two X-ray pulses vanishes around 10 keV 
(confirmed with RXTE-PCA) and a new pulse appears above 
20 keV {Sect. \3~2J\ and Fig. 0. 

3) Three XMM-Newton observations in 2003-2004 give sta- 
tistically identical pulse profiles, but the profile measured by 
ASCA-G IS in 1999 following a possible glitch of 4U 0142+61 
(see also iMorii et al.l l2005bt iDib et al.1 120071) is significantly 
different with different relative strengths and spectra (0.8- 
10 keV) of the two pulses visible in the pulse profiles below 
10 keV (Sect. \TZ2\ and Fig. 0&©. 



4.4. Pulsed spectra 

1) ISGRI measures above 20 keV a total-pulsed spectrum up to 
160 keV with a shape consistent with a power law with index 
0.40±0.15, significantly different from that of the total spec- 
trum. This implies that the difference spectrum, the DC com- 
ponent, is a genuinely different component (Sect. 13.3.71 Fig. \8\ 
and Tabled . 



2) Below 10 keV, there is evidence for long-term time vari- 
ability. The total pulsed fluxes measured for the 0.8-2.0 keV 
band with ASCA-GIS in 1999 and in the three XMM-Newton 
observations in 2003-2004 are comparable, but in the 2.0- 
10.0 keV band different flux values are measured. The ASCA- 
GIS value is significantly lower by ~35% than the mutually 
consistent XMM-Newton values, and the RXTE-PCA flux ap- 
pears to be in between. The ~35% decrease in total pulsed flux 
in the ASCA-GIS observation can largely be assigned to a de- 
crease in the strength of only one of the X-ray pulses (Sect. \3.3\ 
FigMand Tablesffl&M. See also lDibetal1(l2007b for a signif- 
icant increased pulsed flux of ~3 6% since 2004 and d uring the 
bursting activities in 2006-2007 dGonzalez et al.ll2008h . 

3) The total pulsed fraction increases from ~10% at 1 keV to 
possibly as high as ~40% in the 40-80 keV band. Above these 
energies this fraction is uncertain, but might become close to 
100%, depending on the actual shapes of the total and pulsed 
spectra (Sect. \3.3.3\ and Fig. \TP\ . 

4.5. Phase-resolved pulsed spectra 

1) Phase-resolved spectroscopy over the total high-energy band 
reveals the identification of at least three genuinely different 
pulse components with different spectra: a) a component in 
phase interval Phlll (phase 0.85-1.35) peaking at phase 1.1, 
present and dominant from 0.8 keV up to 160 keV. b) a hard- 
X-ray component in phase interval Ph I peaking in the ISGRI 
band above 20 keV at phase ~0.55. c) a soft X-ray component 
in phase interval Phi + Phil (phase 0.35-0.85) peaking in the 
profiles for energies below 10 keV at phase ~0.65, and not vis- 
ible above 10 keV. d) the soft X-ray component mentioned in 
c), appears to consist of two parts, a softer contribution from 
phase interval Ph II than from phase interval Ph I, consistently 
seen for all reported X-ray observations (Sect. \3~~3~\ Fig.^\and 
Tables[7\&^. 

2) The total pulsed spectra and that of Ph III measured with 
RXTE (1996-2003) and consecutively with ISGRI (2003- 
2007) connect smoothly, and the phase-resolved spectra mea- 
sured with ASCA (1999), XMM-Newton (2003-2004) and 
RXTE (1996-2003) are relatively similar. This is indicative 
of a remarkably stable geometry of the contributing emission 
sites at the surface and/or in the magnetosphere of the magne- 
tar. There is no evidence for a drastic reordering (Sect. 13.3.71 
Fig. \8\and Tabled . 

5. Discussion 

In this section, we discuss how these new findings fit within 
the theoretical models attempting to explain hard X-ray emis- 
sion from AXPs. Before the discovery by INTEGRA L of non- 
therm al emission from AXPs above 10 keV. ICheng & Z hang 
(2001) modeled the production of high-energy gamma radia- 
tion in vacuum gaps proposed to be present in the outer magne- 
tospheres of AXPs. Gamma-ray emission at the polar caps will 
be quenched due to the strong field, but far away from the pul- 
sar surface gamma radiation could be emitted because the local 
field will drop below the critical quantum limit. This radio- 
pulsar scenario could offer a stable configuration, suggested 
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by our results to be required, but their model calculations do 
not reproduce the hard spectra and high X-ray luminosities 
found for e.g. 4U 0142+6 1 . After the discovery by INTEGRAL 
three approaches were elaborated specifically to address the 
non-thermal, very luminous hard X-ray emi ssion from AXPs: 
1) a quantum electrodynamics model by Hevl & Hernquist 
(2005a b); 2) a corona model by Beloborodov & Thompson 
(200jand most recently bv iLvubarskv & Eichlerl J2008I): 3 ) 
a resonant upscattering model by tearing & Harding! (12007b . 
Below, we will consider each of these attempts. 



5.1. Quantum electrodynamics model 



Hevl & Hernquisll (I2005al) created a model to explain the 
origin of the SGR and AXP bursts through fast Quantum 
Electrodynamics (QED) wave modes which could be triggered 
when a starquake or other rearrangement in the magnetic field 
occurs. In their fast-mode break down model, magnetohydro- 
dynamic (MHD) waves generated near the magnetar surface 
propagating outwards through the magnetosphere are altered 
by vacuum polarization. Hydrodynamic shocks can be devel- 
oped and electron-positron pairs produced. An optically thick 
fireball could be created emitting mainly X and gamma rays, 
which can be observed as a SGR or an AXP burst . To ex- 
plain stable non-thermal emission iHevl & Hernquisll (|2005b) 
extended this model to weaker fast modes that do not cre- 
ate an optically thick fireball. The resulting spectrum depends 
on the total energy delivered by the fast mode and the maxi- 
mum magnetic field strength. The latter value determines two 
energy values where the spectrum breaks. The peak energy 
in the spectrum is given by Ebreak ~ 1600-^mc 2 , where 
Bait is the quantum critical field of 4.4 x 10° G, and B the 
neutron-star magnetic field. The other break energy is given 
by Eo ~ 10~ 6 Y~mc 2 . The spectral shapes below, between and 
above these break energies are power-law like with photon in- 
dices of 0, 1.5 and 3, respectively. The power law below Eq 
can extend down to the optical regime, in which pulsed non- 
thermal optical emission (Eichler et al. 2002) has been detected 
from4U 0142+61 dHulleman et alJ2004tlKern & Martinl2002t 



Dhillon et al.ll2005h . 



Peak energies predicted by this model all exceed 1 MeV. 
The fast mode cascade of pairs is not effective if the energy in 
the pairs (£break) is less than 2mc 2 . This natural limit is called 
the minimal model where £break ~ 1 MeV. In this case, the 
pairs are created in a (local) magnetic field much weaker than 
flrrit. lHevll (12007) examined the angular dependence of the non- 
thermal emission in the weak field regime (Bns ^ Sent)- The 
produced spectrum strongly depends on the angle between the 
magnetic field and the wave propagation direction. Specifically, 
the non-thermal hard X-ray emission is expected to be largest 
when the observer's viewing direction to the pulsar is perpen- 
dicular to the local magnetic field, while the thermal emission 
is expected to peak when the observer views the polar regions. 
The difference between the phases at which the thermal and 
non-thermal emissions are peaking can range between 90° and 
180° if the angles a + B > 90°, where a is the angle between 
the rotation axis and the observer's line of sight and B the an- 



gle between the rotation axis and the magnetic (dipole) axis. If 
a + B < 90°, then the phase difference between the non-thermal 
and the thermal emission is 180°. 

Comparing the above QED model predictions with our re- 
sults, we note that the power-law spectral index of the total 
emission above 20 keV of 4U 0142+61, 0.93±0.06, is harder 
than the predicted value of 1.5. Furthermore, the estimated 
break energy of ~270 keV is far below the allowed minimum 
value of ~ 1 MeV, and would mean a neutron star magnetic field 
strength of 1.3 x 10 16 G The spectral index of the total pulsed 
emission of 0. 40+0.15 is inconsi stent with the value 1.5. We 
note that Hevl & Hernquistl (2005 a) do not make a distinction 
between the pulsed and total emission. 

Finally, there is a clear prediction of a phase difference 
between the thermal and non-thermal emission: Our phase- 
resolved spectroscopy shows that the actual situation is much 
more complex: one pulse exhibits a spectrum composed of a 
thermal and a non-thermal component in the same phase inter- 
val (Ph III); the thermal pulse in phase interval Ph I+II seems to 
vanish above 10 keV and a non-thermal pulse shows up above 
20 keV in Ph II, also at a negligible phase difference. 

5.2. Corona model 



The corona model developed by iBeloborodov & Thom pson 
( 20 07]) is an exten sion of earli er work by 



Thompson & Beloboro dov | (| 2005l) and iThompson et al. 



(120021) . iThompson et alJ ([2002) extended the magnetar model 
by a twisted magnet osphere, which was introduced by 



Thompson et alJ Q2000J) to explain some observational char- 
acteristics seen in SGR 1900+14. Accordingly, the assumed 
dipolar field gets twisted due to continuous seismic activity 
as a result of the extreme internal magnetic field and a quasi 
stable loop (tube) can be formed with a maximum radius of 
~2R NS . 

Essential in this scenario is an electric field parallel to the 
magnetic field which is generated by self induction. This par- 
allel electric field is capable of accelerating particles and initi- 
ates an electron-positron avalanche. In this way, toroidal field 
(magnetic) energy is converted into particle kinetic energy. The 
energy production in the tube is related to the voltage differ- 
ence (<D e ) between the footpoints of the loop. The dissipation 
rate is given by L diss ~ 10 37 A^(B/10 15 G)(a/i? NS )(e«D e /GeV) 
ergs -1 , where A(/> is the twist angle (in radians) and a is the 
size of the twisted region. <?<D e must be sufficiently high for 
charged particles to be lifted from the surface of the neu- 
tron star (~200 MeV) or to create pairs in the magnetosphere 
(~y ± mc 2 - 0.5y ± MeV) which are accelerated to ~0.1-1 GeV 

If the voltage is too low to create a plasma, the tube regu- 
lates itself by either reducing the current or an electric field is 
generated (dE/dt ~ Vx B - j) until it is high enough to provide 
the plasma. If the voltage is too high, the parallel electric field 
will be screened by polarization of the plasma. Through this 
self-regulating voltage 'constant' pair production is achieved, 
feeding the high-energy emission. The dissipation time scale 
of the loop, i.e. the unwinding of the twist, is of the order 1- 
10 years. It is determined by the voltage difference <J> e and the 
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corona luminosity Ldiss- For a global twist of a dipole this be- 
comes: fdecay = 0.8(/?Ns/c 2 )£diss/ ( l ) e 2 - Therefore, the stronger a 
twist (brighter corona) the longer the live time. 

The downward beam of electrons and positrons through the 
tube will heat a transition layer between the atmosphere and 
the corona up to ~200 keV, radiating bremsstrahlung to cool. 
In this thermal transition layer additional pairs can be created 
after electron ion collisions. The emerging spectrum will have 
a power-law shape with index F ~ 1 below the break energy. 
These bremsstrahlung photons can convert into pairs if the pho- 
ton energy exceeds 2m e c 2 / sin 9. In that case synchrotron emis- 
sion will be produced. Either way, the maximum energy of pho- 
tons escaping the strong magnetic field is ~1 MeV. 



to 100 keV, given the pulsed fraction of 
energies. 



50% at these high 



Beloboro dov & Thomps on ( 2007) also address the produc 



tion of non-thermal emission below the INTEGRAL-ISGRI 
energy window. In the 2-10 keV band resonant upscattering 
of keV photons can explain the soft power-law like X-ray 
tai l (>5 keV) to the black -body com ponent as first described 
by Thompson et al. (2002) (see also Lyutikov & Gavriil 2006; 



Gtiver et alj 12006c Fernandez & Thompson! 120071) . For the 



non-thermal optical and NIR emission four possible emis- 
sion mechanisms are discussed with a curvature-radiation 
origin being the most plausible explanation. See also the 
interpretati on of the NIR , optic al and soft X-ray o bser- 
vations by [Hu lleman et al. (2004J), iMorii et al] d2005ah and 
Durant & van Ke rkwiik (2006c). 

A critical aspect in the corona model is the dissipation time 
scale of the twist. This determines how long a quasi stable 
loop can produce a non-thermal component in the pulse pro- 
file. Once the loop vanishes, the pulse profile should drasti- 
cally change shape. When a next starquake produces a new 
tube, a new non-thermal pulse at a different phase should ap- 
pear in the profile. The unwinding of such quasi stable loops 
feeds the non-thermal emission, and it seems plausible that also 
long-term flux variability occurs. Our results give increasing 
evidence that the non-thermal emission and the global struc- 
ture of the pulse profile are stable over many years. During the 
INTEGRAL observations covering more than three years, we 
do not see significant variation in intensity and spectral shape 
of the total and pulsed emission above 20 keV. Furthermore, 
the RXTE monitoring observations covering the eight years be- 
fore the INTEGRAL observations give time-averaged spectra 
for the total and total pulsed emission, as well as for the spec- 
trum of the main pulse (Ph III) in the pulse profile that connect 
smoothly to the INTEGRAL-ISGRI spectra. In addition, the 
pulse maxima and the minima in the pulse profiles measured 
by ISGRI and RXTE line-up, suggesting a very stable geome- 
try and intensity at hard X-rays over more than 10 years. This 
seems to be at the limit of what the corona model still can sus- 
tain for a bright corona. 

The proposed bremsstrahlun g spectrum with index F ~ 1 
dBeloborodov & Thompsonll2007l) is consistent with our total 
spectrum for 4U 0142+61 wifhr = 0.93 ±0.06 below the break 
energy. However, we showed that the total pulsed spectrum is 
much harder with F = 0.40 + 0.15, and, as a consequence, the 
DC difference spectrum much softer, both not consistent with 
a bremsstrahlung origin. The DC spectrum extends at least up 



Lvubarskv & Eichle r (2008) revisited the heating by the 



downward beam of electrons and positrons of the transition re- 
gion and conclude that this mechanism is incapable of heating 
the atmosphere due to suppression of the required two-stream 
instability. They propose, rather, that a hot, tenuous atmosphere 
forms with a temperatu re of ~ 1-2 MeV, an order of magn itude 
larger than inferred by Thompso n & Beloborodov! (J2005) and 
Beloborodov & Thomps on (2007). In this hot atmosphere hard 
radiation is generated via bremsstrahlung. Furthermore, pairs 
are easily produced and fill the whole magnetosphere forming a 
hot corona. In the latter collisionless interaction of the primary 
beam with the pair plasma in the corona heats the pairs even 
more. Most of this coronal energy is released in the hard X-ray 
and soft gamma-ray bands by Comptonization, so that the over- 
all spectrum of the source is a superposition of bremsstrahlung 
radiation from the hot atmosphere and Comptonization radi- 
ation from the extended corona. The latter radiates in all di- 
rections, and the spectrum extends to the MeV band. The 
Comptonization radiation dominates below ~100 keV with 
photon spectral slope of F ~1— 2, with a bremsstrahlung spec- 
trum (F ~ 1) above this break energy. This overall spectral 
shape is not consistent with our results on 4U 0142+61. The 
Comptonization spectral shape below the break energy might 
be marginally consistent with that of the total spectrum, but 
by far too soft for the total pulsed spectrum. Furthermore, a 
bremsstrahlung spectrum after the break energy would mean 
an upturn of the spectrum, extending up to the MeV band. This 
is inconsistent with the INTEGRAL-SPI results above 100 keV 
and more importantly, with the COMPTEL upper limits in the 
MeV band. Finally, it is hard to imagine how a Comptonization 
component originating in an extended corona radiating in all 
directions can produce the phase distribution of 4U 0142+61 
with well defined pulses and, for 4U 0142+61, a sharp mini- 
mum. 

Our concerns on the long-term stability/variability 
of the corona model a s expr essed for the work by 
Beloborodov& Thompson (200 71), are also valid for the 
approach of iLvubarskv & Eichleri (2008). 



5.3. Compton upscattering model 



Baring & Hardingl (120071) consider resonant, magnetic 
Compton upscattering to explain the hard X-ray emission. 
Ultra-relativistic electrons, accelerated along either open or 
closed magnetic field lines, upscatter the thermal soft X-ray 
atmospheric photons. This process is very effective close 
to the magnetar surface (<10/?ns) where the magnetic field 
is still strong. For high-B isolated pulsars it is shown that 
acce leration of electrons to ultra-relativistic energies can occur 



(e.g.lSturnerll995tlHarding & Muslimovll998tlDvks & Rudak 
20001) . However, if AXPs with their long periods have also 
a dipolar structure at the surface, the acceleration zone cap 
radii would become extremely small (of order 50 m). As a 
result, the relativistic electrons would have to be concentrated 
in such a narrow column with charge densities far exceeding 
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the G oldreich-Julian maximum density (Gol dreich & Ju lian 
19691) . They argue, however, that non-dipolar structure at 
the surface is energetically feasible, allowing the electron 
acceleration zone to cover a much larger range of colatitudes 
than assigned to a standard polar cap. 

The Compton scattering formalism invokes full relativis- 
tic and relativistic QED effects in high magnetic fields. 



Referri ng to earlier studie s by i.e. iDa ughertv & Harding 



dl986l) : iHo & Epstein] d 1989k and iGonthier et al.1 (l200Ci . and 



applications in the magnetic Thomson lim it for old gamma- 
ray burst sce narios (Dermer 1990; Baring 199 4]) and for pul - 
sar contexts JDaughertv & Hardindll989l:ISturner et al.|[l995l) . 



Baring & Harding (2007) make a first attempt to calculate ge- 
ometries and emission spectra produced in the 'resonasphere' 
of AXPs for a simplified scenario. They consider the case for a 
magnetic dipole field geometry and an injected monoenergetic 
electron distribution as well as monoenergetic incident photons 
of energy commensurate with the thermal photon temperatures 
kT ~ 0.5 - 1 keV observed in AXPs. Important first conclu- 
sions are that the Compton resonasphere for long period AXPs 
is confined to within few stellar radii of the surface at higher 
field locales, and that the emission spectra are considerably 
softer than the hard X-ray tails seen in AXPs, and extend to 
energies higher than allowed by the COMPTEL upper limits, 
as e.g. shown in Fig. [3] The spectra do break above a maxi- 
mum photon energy which is dependent on the magnetic field 
and the Lorentz factor. Furthermore, the observed spectra de- 
pend very sensitively on the viewing direction of the observer 
with respect to the magnetospheric geometry. For investigat- 
ing such predictions, phase-resolved spectroscopy as applied in 
this work offers a very sensitive diagnostic tool. However, more 
complicated modelling addressing non-dipolar field topologies 
and more realistic electron and photon spectra are required be- 
fore our detailed findings can be confronted with the predic- 
tions for hard X-ray emission from resonant Compton upscat- 
tering in AXPs. What ca n be noted, is that the pre liminary ex- 
ample spectra derived bv lBaring & Hardind ( 12007b are too soft 
and extend to energies much higher than can be permitted by 
the COMPTEL upper limits. 



5.4. Concluding remark 

The above discussions indicate that our detailed results as sum- 
marized in the previous section require, for all attempts to 
model the non-thermal emission above 20 keV, consideration 
of a full three-dimensional geometry taking the relevant angles 
(magnetic axis with respect to spin axis, viewing angle) into 
account, as well as the physical production processes taking 
place on the surface of the neutron star and in different sites 
of the atmosphere and magnetosphere. Such analysis would be 
similar to what is done to model the non-thermal emission from 
radio pulsars. Our results seem to indicate that the non-thermal 
persistent emission from AXPs is persistent, requiring a sce- 
nario with a stable geometry. To obtain further information on 
the geometry, including more information for the different rele- 
vant angles, similar detailed results for the other non-thermally 



emitting AXPs (see iKuiper et al J 120061) are important to give 
further constraints on the theoretical modelling. 
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